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Fermentation is generally modelled by kinetic equations giving the time evolutions for biomass, 
substrate, and product concentrations. Although these equations can be solved analytically in simple 
cases if substrate/product inhibition and biomass death are included, they are typically solved 
numerically. We propose an analytical treatment of the kinetic equations — including cell death 
and an arbitrary number of inhibitions — in which constant yield needs not be assumed. Equations 
are solved in phase space, i.e. the biomass concentration is written explicitly as a function of the 
substrate concentration. 

Keywords: fermentation, growth kinetics, model, inhibition, biomass yield 



Introduction 

Several models have been proposed to describe 
the kinetics of fermentation, giving the time evolu- 
tions for microbial mass X , substrate S, and prod- 
uct P. Marin and Mitchell et al. have re- 
cently reviewed them. Kinetic equations are gen- 
erally of the form 



dS 1 dX 

= TO A, 

dt Y dt 



dX _ S 
~dT ~^Y+K~i 
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where to is the maintenance coefficient 0, 0i A* 
is the maximum specific growth rate and K\ is 
Monod's constant 5j. Equation @, originally pro- 
posed by Monod j5j, can be modified to account 
for product inhibition [1 IS 0, H 0, EH , substrate 
inhibition [S ITTl IT^ , and biomass death [l3l IT^| . 
Other models are possible; Cramer et al. for 
instance did not correlate ethanol production to 
biomass growth. 

The biomass yield, —dX/dS, is often assumed 
to be constant [S 0, 0, 0, 0, . In this case 
X - X° = Y (S - S°) where X° and S° are the 
initial values of X and S. Equation J3J) then gives 



/j, t = In 



X 



X° S° + X°/Y \ X 



In- 



X 



hi 



However biomass yield is not necessarily con- 
stant. Thatipamala et al. |l9l | for instance found 
that it decreased from 0.16 to 0.03 when ethanol 
increased from to 107 g/L. In such cases the 
time dependences of 5* and X are not obtained 
analytically, instead Eqs. QJ and J2J) are typically 
solved numerically. However since monitoring of 
the biomass concentration can be used to follow 
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the evolution of fermentations an explicit expres- 
sion linking X and S could be very useful. 



Complete kinetic equation 

Equation @ assumes that the microorganisms 
are not affected by their environment. In actual- 
ity fermentation can be slowed by large amounts 
of substrate (substrate inhibition) or of product 
(product inhibition). Biomass death may also oc- 
cur. In this section we write a more complete ver- 
sion of Eq. (J2J in order to take these effects into 
account. 

Both substrate and product inhibitions can be 
modelled by an equation of the form 



dX 
~dt 



S 



-x. 



(3) 



The {Ki} can be either positive or negative de- 
pending on the nature of the inhibition (K\ is al- 
ways positive, it is Monod's constant). They can 
refer to different kinds of sugars, nutrients (e.g., 
nitrogen) or products. A 4 '/ 11^=1 (1 + S/Ki) must 
always be smaller than 1 since inhibition slows 
growth and it is equal to 1 in the absence of in- 
hibition. // is given by 



n 

ki<0 



where p° = if /(c^S ) and k t = K./S . The prod- 
uct is over all i corresponding to product inhibi- 
tion, ai is the stoichiometric coefficient for product 
Pi. (1 + S/K^il + S/K^) can refer to two different 
substrates or products but it can also refer to a sin- 
gle variable if it is found that (1+S/Ki)(l+S/Kj) 
can provide a better fit to the data than 1 + S/Ki 
alone. 
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Inhibitions can slow fermentation but they can- 
not make the biomass decrease. In order to ac- 
count for a decrease in viable biomass, cell death 
must be included in the model. Equation (J2J) be- 
comes liH 



x/s 



dX 

~dt n 5 + A 



Kr 



-X 



(4) 



where Arj is a constant (g/L h). The first term 
is biomass growth and the second term is biomass 
death. The way Eq. QJ was written dX/dt was 
the same as biomass growth. But in order to 
use Eq. J2J one needs to distinguish between the 
change in the biomass due to growth (which will re- 
sult in sugar consumption) and the effect of death 
(which has no effect on 5). dS/dt depends only 
upon the growth part of dX/dt. 

The complete equation accounting for inhibi- 
tions and cell death is obtained by merging Eqs. © 
and (D 



dX p,p! S 



dt Kx J-L 

2=1 



X- Pf ,J2- x > (5) 



where m = mY Kn{i)/(oii 5°). 
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FIG. 1: x/s as a function of s for ip = 10 2 (solid line), 
ip = 10 _1 (dashed line), and ip = 1 (dotted line). 



x and s are normalized biomass growth and sub- 
strate consumption respectively. They are zero at 
t = and they go to 1 when t — + oo. s depends 
only on x and on ip. Figure ^ shows x/s as a func- 
tion of s for three values of ip. When ip -C 1 (solid 
line) x « s, but the assumption of constant yield 
is less good for larger values of ip. 

A similar result can be obtained from Eqs. 
and JSJ, accounting for product and substrate in- 
hibitions combined with cell death. Equation JJjJl 
becomes 



Solution in phase space 

In this section we will show how one can ob- 
tain X as a function of 5 from Eqs. and (JSJ, 
without assuming constant yield. But first we will 
proceed from Eqs. (JIJ and in order to describe 
the procedure on a simpler example: 

dS _ dS/dt 1 mS + Ki 
dX ~ dX/dt ~ ~Y ~ Y 5 

which gives 



1 + pdX _ S°ip 



Y dS ' S + S°ip' 
The dimensionless constants p, ki, and ip are 



s° 



1+p 



Letting 



YS° 



s° 



we obtain 



a; = s + ip In 1 — 



1 + tp 



(6) 



p=^—, h = -^, and ip = P fci. (7) 



a; = TTT^n (1 + p) and s = on , (8) 



(9) 



1 dX 



s+sy^ 



y d5 s , pKi 



/'' 



A", 



, (io) 



where Pi = Pi/faS ) = y\ + 1 - 5/5°. Whereas 
Eq. © was easily integrated, Eq. (|10|l is trouble- 
some. The fraction must be transformed into a 
sum of terms of the form 1/(5 — 5*). 

Theorem 1. Let n G IV, n > 2. Le£ {crj e K" . 
Le£ Q a polynom of degree strictly smaller than n. 



Q(s) 



Q((Ti) 



= E- 

U(s-a,) «=i (s-^n^-^) 



(11) 



Proo/. Let Q = {j e IN*, j < n}. Let {R l } ie n 
the polynoms defined over 1R by Ri(s) = 

U&i (« ~ oj-)- Let (*. *) e Rii^/RiWi) = 1 
if i = I and otherwise. Let T the polynom defined 
over 2R by T(s) = Q(s)-£" =1 Q{°i)Ri{s)l '-Rifo). 
V^? <E fi, T(ae) = 0. T is of degree smaller than or 
equal to n — 1 and it has n roots. It is the zero 
polynom. Hence Eq. (|llfl . □ 

Letting {ci} the roots of /i'(l — s) + 
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pfcillfcil 1 + (1 _ s )l^i\-> we have 



x =// 2J (o» - 1) <t>i In 



™ D 1 / c 2 



(12) 



where 



rl D i 

(Ti + Pi 



(13) 



One should note that the {ci} do not depend on 
cell death. 

In Eq. I|10[l the numerator of the fraction is of 
degree 2 (of degree 1 in the absence of cell death) , 
the denominator is of degree n, with n > 1. Equa- 
tion l|12l) applies only if the degree of the numer- 
ator is smaller than that of the denominator (the- 
orem 1). In the absence of cell death, Eq. I|12fl 
holds if there is at least one inhibition and two 
inhibitions or more are necessary if cell death is 
taken into account. 

Let us consider one inhibition (parameter k 2 ) 
and no cell death. If k 2 large in absolute value, 
S\ ~ kik 2 /ip, s 2 ~l + V>, and fa ~ (— 1)\ 
Thus (pi — 1) 0i In 1 1 — s/si\ ~ s and one recov- 
ers Eq. ©. It is not possible to write x — 
fi'(si — 1) fa In 1 1 — s/si\ because Eq. (fT2*|l applies 
only for at least one inhibition. 
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FIG. 2: Normalized biomass yield, x/s, as a function of 
s, accounting for inhibition, ki = —1.5 (dashed line), 
k'2 — —1.05 (dotted line), k^ = 0.05 (dot-dashed line), 
and ki — > ±oo (i.e. no inhibition, solid line). 
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FIG. 3: Normalized biomass yield, x/s, as a function 
of s, accounting for cell death, k is 0.5 (dashed line), 
0.1 (dotted line), and 0.05 (dot-dashed line). The solid 
line corresponds to the absence of cell death (ft — > oo). 



Applications 



Equation i|12|l is very general as it can account 
for an arbitrary number of both substrate and 
product inhibitions as well as biomass death. In 
this section we consider particular cases: inhibi- 
tion without death, cell death without inhibition 
before combining the two effects. 

Figure |21 shows x/s in the case of one inhibition 
(parameter A^). As expected substrate inhibition 
(dot-dashed line) affects the fermentation from the 
beginning whereas product inhibition matters at 
the end of fermentation. The total yields for 
&2 = 0.05 (dot-dashed line) and k 2 = —1.05 (dot- 
ted line) are similar although the routes are quite 
different. k 2 cannot be between —1 and because 
of its definition, k 2 = — 1 corresponds to K' = if 
one writes (1 - P/K')- 1 instead of (1 - S/K 2 )-\ 
The values of the constants used in Figs. [21 and 
subsequent are [2(|: ju = 0.11 h _1 , m = 0.01 h _1 , 
Y = 0.1, and K = 112 g/L. The initial substrate 
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FIG. 4: Normalized biomass yield, x/s, as a function 
of s. The solid line corresponds to the absence of in- 
hibition and cell death. The dashed line is for product 
inhibition alone (fe = —1.5) and the dotted line for 
biomass death alone (k = 0.1). The dot-dashed line 
shows the combination of the two effects. 
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concentration is S° = 100 g/L and a = 1/2. This 
gives tp k 1 %. 

FigureOlshows x/s as a function of s if cell death 
occurs but there is no inhibition. If n is small, 
s may not reach 1 (stuck fermentation) as shown 
by the dot-dashed line in Fig. [31 Cramer et al. 
|15| found a value for i^D close to 5 kg/L h which 
corresponds to k ks 0.1. 

Figure0]shows x/s as a function of s when both 
product inhibition (parameter kz) and cell death 
(parameter k) occur. The initial product concen- 
trations, {p®}, are all set to zero. The effects of 
product inhibition and cell death do not simply 
add up. Product inhibition alone (dashed line) re- 
duces the yield only slightly. However the yield 
decreases significantly when product inhibition ex- 
ists on top of cell death (dot-dashed line compared 
to dotted line). The same final value of x can be 
obtained with k = 0.1 and hi = —1.5 or with 
K w 0.049 and no inhibition. The dot-dashed line 
in Fig. |U (k = 0.1 and ki = —1.5) is indeed very 
similar to the dot-dashed line in Fig. El (k = 0.05, 
no inhibition). 

Conclusion 

By eliminating time from the kinetic equations, 
one can solve them in phase space, i.e. the biomass 
concentration is obtained as a function of the sub- 
strate concentration. It is then possible to obtain 
an analytical expression for the effect of an arbi- 
trary number of substrate and product inhibitions 



combined with biomass death, without assuming 
constant yield. Once parameters are obtained for a 
certain process, substrate and product concentra- 
tions can be obtained analytically from the mea- 
sured biomass concentration in order to follow the 
evolution of the fermentation. 



Nomenclature 



K, 


inhibition coefficient 


(g/L) 


kj 


dedimensionalized Ki, ki = Ki/S° 


(-) 


K u 


cell death coefficient 


(g/Lh) 


m 


maintenance coefficient 


(h- 1 ) 


"D 


number of biomass death terms 


(-) 


m 


number of inhibitions plus one 


(-) 


Pi 


product concentration 


(g/L) 


Pi 


dedimensionalized Pi, pi = Pi/(ctiS°) 


(-) 


s 


substrate concentration 


(g/L) 


s 


normalized substrate, s = 1 — S/S° 


(-) 


t 


time 


00 


X 


viable biomass concentration 


(g/L) 


X 


normalized biomass, x = (X- X°)(l+p)/(Y S°) (-) 


Y 


biomass yield 


(-) 


o-i 


stoichiometric coefficient for product Pi 


(-) 


K 


dedimensionalized Kr>, k = mY Kr> / (aS°) 


(-) 


A< 


maximum specific growth rate 


(h-i) 


v! 


change in /i due to inhibition 


(-) 


P 


ratio of time constants, p = m Y/ p, 


(-) 


v, 


a root of fi'(l - s) + pfcilljl 1 + (1 — s )/ k j] 


(-) 


<Pi 


defined in Eq. Ilijl 


(-) 


V* 


distance to constant yield, = kip/(X + p) 


(-) 



Superscripts 

initial value 



[1] R. Marin, Am. J. Enol. Vitic. 50, 166 (1999). 
[2] D. A. Mitchell, O. F. von Meien, N. Krieger, 

and F. D. H. Dalsenter, Biochem. Eng. J. 17, 15 

(2004). 

[3] F. Acevedo, CRC Crit. Rev. Biotechnol. 6, 309 
(1987). 

[4] G. A. Hill and C. W. Robinson, Chem. Eng. J. 44, 
B69 (1990). 

[5] J. Monod, Recherche sur la croissance des cultures 

bacteriennes (Hermann, 1941). 
[6] C. D. Bazua and C. R. Wilke, Biotechnol. Bioeng. 

Symp. 7, 105 (1977). 
[7] I. Holzberg, R. K.Finn, and K. H. Steinkraus, 

Biotechnol. Bioeng. 9, 413 (1967). 
[8] S. Aiba, M. Shoda, and M. Magatani, Biotechnol. 

Bioeng. 10, 845 (1968). 
[9] T. K. Ghose and R. D. Tyagi, Biotechnol. Bioeng. 

21, 1401 (1979). 
[10] G. K. Hoppe and G. S. Hansford, Biotechnol. Lett. 

4, 39 (1982). 

[11] C. S. Ough and M. A. Amerine, Am. J. Enol. Vitic. 



14, 194 (1963). 

[12] G. P. Casey and W. M. Ingledew, CRC Crit. Rev. 
Biotechnol. 13, 219 (1986). 

[13] J. E. Bailey and D. F. Ollis, Biochemical engineer- 
ing fundamentals (McGraw-Hill, New York, 1977). 

[14] I. Caro, L. Perez, and D. Cantaro, Biotechnol. Bio- 
eng. 38, 742 (1991). 

[15] A. C. Cramer, S. Vlassides, and D. E. Block, 
Biotechnol. Bioeng. 77, 49 (2001). 

[16] G. R. Cysewski and C. R. Wilke, Biotechnol. Bio- 
eng. 18, 1297 (1976). 

[17] G. R. Cysewski and C. R. Wilke, Biotechnol. Bio- 
eng. 19, 1125 (1977). 

[18] B. L. Maiorella, H. W. Blanch, and C. R. Wilke, 
Biotechnol. Bioeng. 25, 103 (1983). 

[19] R. Thatipamala, S. Rohani, and G. A. Hill, 
Biotechnol. Bioeng. 40, 289 (1992). 

[20] H. H. Topiwala and C. G. Sinclair, Biotechnol. 
Bioeng. 13, 795 (1971). 



